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Abstract 

Usually it is believed that the Hartree-Fock state of quark matter is a Fermi 

gas state with no polarisation of spins. We examine the possibility of the 
>. 
^^ ' polarised quark liquid interacting with the one-gluon-exchange interaction. 

f^ ' It is suggested that the Hartree-Fock state shows a spontaneous magnetic 

^H ! 

On ' instability at low densities through the same mechanism as the appearance of 

ON 

r~| ' ferromagnetism in electron gas. Metastability of the polarised quark liquid is 

Q^' also discussed. 

^ ■ PACS: 04.40.Dg, 07.55.Db, 12.38.Bx, 26.60.-Fc, 97.60.Jd 
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The physical origin of the strong magnetic field in neutron stars is a long-standing prob- 
lem |I|. Recent discovery of "magnetars" [^Q - neutron stars with magnetic field strengths 
in excess of Bcr = mlc^/eh ~ 4.4 x lO^^G and typically ~ lO^^G - also stimulate our interest 
in the origin of their strong magnetic field. It has been widely believed that the magnetic 
field of neutron stars is inherited from the progenitor stars |I| . However, the extraordinary 
magnetic field in magnetars seems to enforce our reconsideration of its origin. Inside neutron 
stars, hadronic matter exists beyond the nuclear density (no — 0.16fm~^), so that it may be 
interesting to consider the hadronic origin of the magnetic field. 

In this Letter we examine the possibility of the spin-polarised quark liquid by the use 
of the Fermi liquid theory, where quarks interact on each other by the one-gluon-exchange 
(OGE) interaction. Since the system should be totally color singlet (neutral), there is left 
only the exchange energy within the Hartree-Fock approximation. It has been suggested 
that three-fiavor quark matter (strange quark matter) can be stable or metastable for a 
reasonable range of QCD parameters |^|-0]. The appearance of strange (s) quarks is favored 
at and above the nuclear density due to the reduction of the kinetic energy of the system and 
the specific feature of OGE interaction; the exchange energy is negative for massive strange 
quarks at low densities, while it is always positive for massless u and d quarks. This feature 
is very similar to the Coulomb interaction for electron gas. In the following we consider the 
possibility of ferromagnetism in strange quark matter on the analogy with electron gas. 

In 1929 Bloch first pointed out the possibility of ferromagnetism of electron gas at low 
density [Q: electrons interact on each other by the exchange Coulomb force, which con- 
tribution is attractive and becomes dominant over the kinetic energy at low density due 
to their different density dependence. Then the unpolarised electron gas suddenly turns 
into the completely polarised state at the critical density (the first order phase transition). 
The critical density is determined by the mass of electrons and the QED coupling constant 
a = e^/47r; lower density, where the nonrelativistic approximation works very well, is favored 
for the spontaneous magnetisation. This idea has been confirmed theoretically by the use 
of the Green function Monte-Calro method M. 



Therefore if the exchange energy is attractive and dominant over the kinetic energy 
at some density, there should appear ferromagnetism also in strange quark matter, as in 
electron gas. The idea can be put into a relativistic formulation utilizing the polarisation 



density matrix |T0[. Then we shall see a difference between theories of electron gas and 
quark liquid. The scattering amplitude of two quarks becomes momentum dependent, and 
includes not only the spin-nonflip term but also the spin-flip one due to the lower component 
of the Dirac spinor, the latter of which vanishes in the non-relativistic limit. Also, OGE 
interaction itself include the spin-dependent (magnetic) interaction as well, which becomes 
effective in the relativistic regime, besides the Coulomb-like (electric) interaction. These 
features somewhat change the mechanism of polarisation in quark liquid from in electron 
gas. Actually we shall see that there is a possibility of spontaneous magnetisation even in 
the relativistic regime, where the exchange energy contribution is no longer attractive, by a 
different reason from electron gas. 

Consider the spin-polarised quark liquid with the total number density of quarks Uq [] ; 
we denote the number densities of quarks with spin up and down by n+ and n_, respectively 
, and introduce the polarisation parameter p by the equations, 

under the condition < p < 1. We assume as usual that three color states are occupied to 
be neutral for each momentum and spin state. The Fermi momentum in the spin-polarised 
quark matter then is 

kt = kF{l + pf'\ kp = kpil - pf'^ (2) 

with the Fermi momentum for the unpolarised matter, kp = (yr^rig)^/'^. The kinetic energy 
density is given by the standard formula. 



^ We, hereafter, consider one flavor quark matter, since OGE interaction never mix flavors. 

^The meaning of spin up and down should be carefully defined in the relativistic theories. Here 
we call them in the rest frame of each particle. Their precise meaning is given later. 
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with the Fermi energy Ep = [m? + kpY^"^. 

The interaction energy is given by the OGE interaction. There is no direct energy 
contribution (Hartree term) and left only the exchange energy contribution (Fock term). 
To characterize the degeneracy of the plane wave solution u'^^\k) [a = 1,2) for a positive 
energy state, we introduce a space-like 4-pseudovector a^^ which in the rest frame is twice the 
three-dimensional mean spin vector (^ with \^\ = 1, a^ = (0,(^). The 4-vector in any frame, 
where the particle is moving with the momentum k, is found by a Lorentz transformation 
from the rest frame, and 

k(C-k) k-C ,,, 

mg{Ek + mq) irig 

with Ek = [nil + Ikp)-*^/^. By choosing C, along the z axis, C, = (0,0, ±1), we see that 



the spinors u^°''{k) are eigenstates of the operator —W ■ a/niq, where the 4-vector W^ i 



is 



the Pauli-Lubanski vector, W^ = —\t^upak^<^^'^- For the standard representation of u^"\k) 
pT| , the eigenvalue is +l/2(spin-up) for u^^\k), and —1/2 (spin-down) for u^'^\k). Then the 



polarisation density matrix p is given by the expression, 

p{kX) = ^i^ + mq)Pia), P(a) = 1(1 + 75^), (5) 

which is normalized by the condition, trp(A;, () = 1. 

Let us consider the interaction between two quarks with momenta, k and q, and the 
different spin vectors, C ^-iid (^', respectively. They are specified by the polarisation density 
matrices, p{kX) and p{q,C)- The Landau Fermi-hquid interaction /kCqC is related to the 
two-particle forward scattering amplitude |I2|, 



Ac,.c' = ^^-^K,.C' (6) 

where A^kCqC is the usual Lorentz invariant matrix element. The color symmetric matrix 
element is given by the exchange term only; the direct term vanishes because the color 
symmetric combinations (~ trAa) does not couple to gluons. Thus 



ML^C = -/^tr(A,/2A,/2)t.(^')(q)7^w(^)(k)«(^)(k)7^«(^')(q: 



= l/itr[7,p(fc,C)7Mg,C')](^^- (7) 

Substituting Eqs. (|) and © into Eq. (0), we finally find 

Mlr^,, = Ig^ \lml - A: ■ g - (k • C)(q ■ O + ^oC • C 
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(^_&fc +mg)[J^g + nig) 

+ (k-q)(C-k)(C'-q)}]^^^. (8) 

after some manipulation of 7-matrix algebra. If we choose both (^ and C' along the z axis 
in parallel, C = C = (0,0,±1), we have the spin-nonfiip amplitude A^^'^^^ ^^, while if we 
choose them in anti-parallel, C = ~C' = (0, 0, ±1), we have the spin-flip amphtude A^kq *^- 
Each form of the spin-nonflip or spin-flip amplitude is complicated, but the average of them 
gives a simple form, 

s 2 „ 2m^ - k ■ q 



which is nothing but the matrix element for the unpolarised case |jT2|. In the nonrelativistic 
limit, mq ^ |k|, |q|, the matrix element is reduced into the form. 



_ 2 , m^(l + C- CO 



•^kcqc = 7;9' '), V > (10) 



so that there is no correlation between quarks with different spins. On the other hand, there 
is some correlation included even in the Hartree-Fock approximation for the relativistic case. 
After summing up over the color degree of freedom and performing the integrals 
of the Fermi-liquid interactions, f\^^°^ "^ = {nig/ Ek){mq/ Eq)Ai1l^^^ ^^ and /^q ^ = 
{nig/ Ek){niq/ Eg)Ai\^j^ *^, over the Fermi seas, we have the exchange energy density e^x con- 
sisting of two contributions, 

_ nonflip I flip /-i -i \ 

tex — tex ' ^ex ■ \^^J 

The spin-nonfiip contribution can be written as 



^nonflip 



Oic 

4vr5 



ilikt) + lik^)), 



where ac is the fine structure constant of QCD, ac = g'^/4:Tc. The integral I{kp){i 
given by 



/(4) = yf. (?. - -"' 



Ki + i{k^p) 



(12) 
±) is 

(13) 



with ^i = Epkp — m? \n{Ep + kp/nig) and Ki = kp — rriq ln{Ep + kp/rriq), and the remaining 



integral 



F ^ {Eg + mg){Eg + 2mg)+ml 



K,{q,k^p) 



qdq 



Eg + m, 



-^2(g,A:» 



with the functions, 



K,{q,k'p) = iE],- Eg) In 



ml - E'pEg + k'pq 



w? - E'pEg - k'-pq 






m„ 



and 



K2{q.K 



ik'^ - q') In 



ml - E'pEg + k'pq 



ml - E'pEg - k'pq 



+ qEg In — ^ -^ 



-2k'pq 



Similarly the spin-fiip contribution can be written as 



ar 



^ex ~ '^TZE'^y^F^kFj- 



The integral J{kp,kp) is given by 



J y Kp , Kp J 



TT 



e+e- 



2m,„ 



{^+K.+^_K+) 



J[kp, kp), 



where 



J y Kp , Itp J 



n'^m^n 



qdq 



Ea + m, 



iEgK,{q,kt) - K^iq^kt)) 



pdp 



iEpK,{p,kp)-K2{p,kp)) 



(14) 



(15) 



(16) 



(17) 



(18) 



/o Ep + mg 
For the unpolarised (p = 0) matter the sum of I and J in Eqs. (|T3|) and (p!8| ) gives 



(19) 
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with the Fermi energy Ep for the unpolarised case, Ep = {m? + ky)^'"^ |[l^JT3|] 



In the nonrelativistic case, ky -C rriq, the kinetic energy density is written as 

The spin-flip contribution becomes tiny in this case, and the dominant contribution for the 
OGE energy density in Eq.(]Tl|) comes from the spin-nonflip contribution, especially / given 
in Eq. (|T^), / ~ — 27r^A;^. Thus we find an attractive contribution, 

^e.--|^{(i+pr/^+(i-pr/1- (22) 

The exchange energy contribution is attractive and exactly given by the Coulomb-like (spin- 
nonflip) interaction of gluons. There is no spin-flip contribution in the nonrelativistic limit. 
The form of the energy densities (^) and (|22|) is the same as in electron gas. It is the differ- 
ence of density dependence between these contributions given in Eqs. ( ^1]) and (|22D which 
may cause a ferromagnetic instability: as density becomes low enough, all the particles sud- 
denly align their spins and the completely polarised {p = 1) matter is energetically favorable. 
This mechanism is firstly discovered by Bloch for electron gas 0; in the charge neutral sys- 
tem, electrons tend to avoid the Coulomb repulsion from each other. If all electrons have 
the same spin state, they can maximally avoid the repulsion by the Pauli principle. On the 
other hand, the kinetic energy of electrons obviously becomes minimum for the unpolarised 
{p = 0) case. We can see that there is a trade-off between two contributions. Thus the 
critical density is determined by the mass and the coupling constant. 

In the relativistic case there are some different features from the nonrelativistic case. 
First, there is a spin- flip contribution due to the lower component of the Dirac spinor even 
for the Coulomb-like interaction. Secondly, the transverse (magnetic) gluons becomes im- 
portant, where the spin-flip effect is prominent. Finally, the density dependence of kinetic 
energy as well as the exchange energy is very different [|14||. Before discussing the general 



case, we consider the relativistic limit, k'^p^ niq] the kinetic energy density behaves like 



^kin 



6k p 

8^ 



{{i + vf" + {i-vY"}. 



(23) 



In this limit, I{kp) ~ ti"^ /2k'^p and J{kp, kp) ~ -n"^ /2k~p^kp^. Thus the exchange energy looks 
like 



a. 



87r= 



k%[{l+pY'^ + {1-pY'^ + 2{l~p 



^2\2/3 



}■ 



(24) 



The kinetic energy is simply an increasing function with respect to the polarisation parameter 
p and takes a minimum at p = 0, whereas the exchange energy is the decreasing function 
and takes a minimum at p = 1. This is due to the characteristic feature of the spin-flip and 
the spin-nonfiip interactions. Both give a repulsive contribution in the relativistic limit. In 
the polarised state (p = 1) there is no spin- flip interaction and the spin-nonfiip contribution 
remains unchanged in spite of the increase of the Fermi momentum. The difference of the 
total energy, ttot = ^kin + ^pot, at p = and p = 1 reads 

4 - 2^/3 ar 



5e = etotip = 1) - etot(p = 0) 






24/3 _ 1 



TT 



(25) 



Hence if ac > 3.23 the polarised state is energetically favorable. Thus, ferromagnetism in 
the relativistic limit has a defferent origin from that in the nonrelativistic case: the spin-flip 
contribution favors the polarised state. 

In Fig.l a typical shape of the energy density is depicted as a function of the polarisation 
parameter p, e.g. for the set ,mg = 300MeV of the s quark and ac = 2.2 as in the MIT bag 
model |T5[ 0. 



Fie.l 



For small p -C 1, the energy density behaves like 



(^tot - (^tot{p = 0) = X P +0{p 



(26) 



^The difficulties to determine the values of these parameters have been discussed in ref. ||6[, and 
we must allow some range for them. 
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with x~^ = XkL + Xpot- X is proportional to the magnetic susceptibihty, and its sign change 
suggests a ferromagnetic transition. It consists of two contributions: the kinetic energy gives 
XkL — kp/(37r'^Ep) (the Pauh paramagnetism), which changes from XtL ~ 0{k%) at low 
densities to Xkin ~ 0{kp) at high densities. On the other hand, the exchange energy gives 



Xpot 



I _ 2ack% 



97r3 



6^^+3W-^^lnf^^±^n|^ + 2t|l+ '"^ 



El ' \ ' ' " \ m, Jj El ' Ell ' ^{Ep + m,) 

(27) 



In the relativistic limit, pp ^ ^q, Xpot behaves like 

Xpot- ^^,kp ^^^kp- ^^^kp [AH) 

from Eq. (p^ , where the first term stems from the spin-nonflip contribution, while the second 
term from the spin-flip contribution. Then we can see that the effect of the spin-flip contri- 
bution overwhelms the one of the spin-nonflip contribution. The potential contribution Xpot 
is always negative , and dominant over Xkln ^^ ^^^ densities, while the kinetic contribution 
Xkin is always positive. If Oc > 37r/2 = 4.7, x becomes negative over all densities. 

For a given set of rriq and Oc, X changes its sign at the critical density, denoted by rid, 
and it clearly shows that quark liquid is ferromagnetic below that density. This is equivalent 
with the sign change of the curvature of the free energy at the origin as a function of the 
polarisation parameter p. Note that the ferromagnetic transition in our case is of the first 
order, so that it is not sufficient to only see the susceptibility. Even above that density the 
ferromagnetic phase may be possible. Actually there is a range, rid < nq < nc2, which is 
specified by the condition s.t. x > ^ind e < 0. Above the density nc2 there is no longer the 
stable ferromagnetic phase. However, the metastahle state is still possible up to the density 
nc3, which is specified by the condition s.t. rj = detot/dp |p=i < 0. In Fig. 2 we depict the 
quantities x"^, ^^ and rj as functions of density, e.g. for the set ,mg = SOOMeV and ac = 2.2. 
The crossing points with the horizontal axis indicate the critical deisities nd,nc2 and rics, 
respectively. 

Fig.2 



We can see that the ferromagnetic instabihty occurs at low densities, while the metastable 
state can exist up to rather high densities. 

Finally we show the critical lines satisfying x~^ = 0,6e = and r/ = in the QCD 
parameter [ac and rriq) plane, which seperate the three characteristic regions for a given 
density. In Fig. 3 they are presented at a density Ug = 0.3fm~^. All the lines have the 
maxima at the intermediate quark mass, and the mechanism of ferromagnetism is different 
for each side of the maximum, as already discussed. If we take ruq = SOOMeV for the s 
quark or mq ~ OMeV for the u or d quark, and ac = 2.2 as in the MIT bag model again 
jTEl , the quark liquid can be ferromagnetic as a metastable state. 



Fig.3 



In this Letter we have pointed out a possibility of ferromagnetism of quark liquid: quark 
liquid potentially bears a nature to be ferromagnetic due to the OGE interaction. We 
have calculated the energy density of the spin polarised quark liquid interacting with the 
OGE interaction in a relativistic framework, and found that the spin- flip and spin-nonflip 
contributions exhibit an interesting feature as a function of the polarisation parameter p or 
density. The ferromagnetic transition in quark liquid is of the first order as in electron gas 
and there is a metastable state even above the critical density. 

We have found that ferromagnetic instability is feasible not only in massive quark system 
but also in light quark system. The mechanism is different between two systems: the similar 
mechanism to electron gas can be applied to the former system, while a novel mechanism 
due to the large spin-flip contribution, which arises as a specific feature of the relativistic 
system, may be applied to the latter system. 

We have seen that the ferromagnetic phase is realized at low densities and the metastable 
state is plausible up to rather high densities for a reasonable range of the QCD parameters. 
Our calculation is based on the Hartree-Fock approximation. So we need to examine the 
higher-order gluon-exchange contributions, especially the ring diagrams, which has been 
known to be important in the calculation of the susceptibility of electron gas ||9|Jl9|l. The 
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ring diagrams have been calculated for the unpolarised quark matter |T^ , [T8[| , but there is no 
calculation for the polarised matter. 

The lifetime of the metastable state is another remaining problem to be studied. It may 
be caused by the highly nonperturbative quantum effects. So it should be proportional to 
the exponential factor at least. Hence it must be difficult to estimate the lifetime reliably. 

If a ferromagnetic quark liquid exists stably or metastably around or above nuclear 
density, it has some implications on the properties of strange quark stars and strange quark 
nuggets 0. They should be magnetized in a macroscopic scale. So strange quark nuggets 
may be no longer a candidate of dark matter in the universe, since they can emit the 
electromagnetic waves like pulsars. 

For a consideration of a possibility to attribute magnetars to strange quark stars in a 
ferromagnetic phase, we roughly estimate the strength of the magnetic field at the surface 
of a strange quark star. The stellar parameters of strange quark stars are similar to those 
for normal neutron stars for the typical mass around Mq = IAMq [|1^,|T^. The magnetic 
dipole moment Mq amounts to Mq = fig ■ {Att/3 ■ r^)nq for the quark sphere with the radius 
Tq, where fiq is the magnetic moment of each quark, /i, ~ fJ^NifJ'N :nuclear magneton~ 
5 X lO^^^erg ■ gauss"^) for massive quarks and lO^/iTv for light quarks. Then the dipolar 
magnetic field at the star surface r = /? ~ lOKm takes a maximal strength at the poles, 

Stt /r \^ 

Bmax = Y ( -| j f^qnq, (29) 

which is order of O(10^^~^^)G for Uq = O(0.1)fm~^, enough for the magnitude required for 
magnetars. 

One may wonder whether such a strong magnetic field modifies the Fermi sea of quark 
liquid itself. We can see, however, that the mean magnetic field working on one quark 
, which is produced by other quarks, vanishes due to the uniform distribution of quarks. 
Hence the self-consistency is retained even in such a case. 
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FIG. 1. Plot of the energy density as a function of the polarisation parameter at Ug = O.lfm"^ 

and Uq = 0.2fm^'^. 
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FIG. 2. Density dependence of x~^: ^^ ^^'^ V- ^^ this case, rid — O.lfm"^, nc2 — 0.15fm~^ and 
Ucs > 0.5fm~^. 
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FIG. 3. Phase diagram in the QCD parameter plane. The region above the sohd curve is in 

a ferromagnetic phase. The region bounded by the sohd and dased curves is in a paramagnetic 

phase, but a metastable ferromagnetic state is possible. The dash-dotted curve shows the boundary 

where the susceptibility is vanished. The dotted curve shows the corresponding one to the solid 

curve with the nonrelativistic formulae given in Eqs. (21) and (22) for comparison. 
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